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The universality of renormalization group limit cycle behavior is illustrated with a simple discrete 
Hamiltonian model. A non-perturbative renormalization group equation for the model is soluble 
analytically at criticality and exhibits one marginal operator (made necessary by the limit cycle) and 
an infinite set of irrelevant operators. Relevant operators are absent. The model exhibits an infinite 
series of bound state energy eigenvalues. This infinite series approaches an exact geometric series 
as the eigenvalues approach zero - also a consequence of the limit cycle. Wegner's eigenvalues for 
irrelevant operators are calculated generically for all choices of parameters in the model. We show 
that Wegner's eigenvalues are independent of location on the limit cycle, in contrast with Wegner's 
operators themselves, which vary depending on their location on the limit cycle. An example is then 
used to illustrate numerically how one can tune the initial Hamiltonian to eliminate the first two 
irrelevant operators. After tuning, the Hamiltonian's bound state eigenvalues converge much more 
quickly than otherwise to an exact geometric series. 
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I. INTRODUCTION 

The renormalization group has been known to be able 
to produce (in principle) a limit cycle for a long time 
A limit cycle is an alternative to a fixed point, al- 
though a limit cycle necessarily implies the existence of 
a fixed point for a discrete version of the renormalization 
group transformation: see below. Recently, the two of us 
have described a discrete, analytically soluble Hamilto- 
nian matrix that can be analyzed using a renormalization 
group procedure. The result, after an initial analysis, was 
a remarkably simple form of limit cycle involving only one 
coupling constant ■ Our model can also lead to chaotic 
renormalization group behavior |2|. Here, we provide a 
more complete analysis of the same model. The more 
complete analysis leads to renormalization group trans- 
formations involving either a finite set of coupling con- 
stants (with 2, 3, 4, or more couplings, as one wishes) 
or in the limit to a transformation in a space of func- 
tions of one variable. The discrete coupling constants 
become identified with the coefficients of the function's 
expansion about the origin of its argument. The focus of 
the more complete analysis is universality and the nature 
of marginal and irrelevant operators (in Wegner's sense 
0) in the model in the presence of the limit cycle. The 
possibility for chaotic behavior is disregarded here. 

The limit cycle picture of the model may be helpful in 
qualitative understanding of more complex systems. In 
particular it is already known that a three-body Hamilto- 
nian with short range forces exhibits a limit cycle when 
there is a two-body bound state at threshold Q|. Sys- 
tems that might be understood approximately through 
the three-body limit cycle example include nuclear [4j, 



and atomic (g three-body systems with short-range inter- 
actions or many-body systems with similar interactions 
that have not been studied yet 0, 0, @- hi the exact 
critical limit of the three-body case there is a three-body 
bound state spectrum that has an infinite sequence of 
states converging to zero energy as a geometric series; 
the infinite sequence of bound states was noticed 0, 
well before an underlying renormalization group limit cy- 
cle was identified. 

Unfortunately, physical three-body systems known to 
date do not have bound states precisely at threshold, 
the requirement for criticality and an infinite sequence 
of bound states. This makes it necessary to understand 
corrections to critical behavior in the three-body case. 
Obtaining these corrections is a very complex undertak- 
ing that is far from complete. One purpose of this paper 
is to provide an analysis of corrections to the limit cycle 
behavior due to irrelevant operators in a much simpler 
model to analyze than is the three-body Hamiltonian al- 
ready under investigation |Tlj . 

A second purpose of our paper is to lay a basis for 
searches across quantum condensed matter physics more 
broadly for evidence of limit cycle behavior. We will 
stress those aspects of limit cycle behavior that are or 
could be universal, and therefore present if examples of 
limit cycles are found in systems other than the three- 
body case already under study. One very distinctive 
phenomenon to be on the lookout for is the presence of 
geometric series of bound state energies (discrete energy 
levels) converging toward zero. One may have to tune 
a real system by applying external fields and stresses 
or introducing impurities in order to approach critical- 
ity. Therefore, our discussion of corrections to criticality 
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leads to an example of how one can tune a finite model 
Hamiltonian to achieve rapid approach of its eigenvalues 
to their limiting behavior for low energies. 

Despite the complexity of a non-perturbative limit cy- 
cle as opposed to a simpler fixed point, the model dis- 
cussed here is largely soluble analytically and even when 
numerical procedures are needed, they demand relatively 
little in the way of computer time or complex computer 
programs. To be more specific, matrix elements of the 



II. 



RENORMALIZATION GROUP 



model Hamiltonian studied here, H n 



2] 



(m\H\n), are 



H. 



gN — ihNSrnn) , 

where m and n are integers. For m = n, 5 mn = 1 and 
Smn — 0. For m ^ n, 5 mn — and s rnn = (m—n)/\m—n\. 
The numbers E n — b n with b > 1, are eigenvalues of 
the operator Hq, which is defined in the same basis to 
have matrix elements (m\Ho\n) — H mn (0,0). The basis 
states are normalized, (m\n) — S rnn . The energies E n 
of the basis states are limited by a cutoff A = b , so 
that m,n < N. The indices can be also limited from 
below by a large negative integer M, in order to make 
the Hamiltonian matrix finite. But the lower bound M 
will be set equal to — oo in most of the discussion that 
follows. 

Hamiltonians defined by Eq. have a general ul- 
traviolet logarithmically divergent structure, with both 
real and imaginary parts contributing to the divergence. 
Therefore, the coupling constants c/jv and hjy are ex- 
pected to depend on the cutoff N = In A/ In b if the phys- 
ical content of the theory, e.g., the energy spectrum, is 
to be independent of the cutoff. In fact, g^ exhibits 
asymptotic freedom as a function of A if hjsi — 0. When 
h-N 7^ OjOat exhibits instead the limit cycle behavior (or 
chaos) ,2). The universality of the limit cycle behavior is 
studied here. 

Section ITTI introduces the renormalization group (RG) 
equation in the model and sets the stage for next sec- 
tions. Section II I II describes the limit cycle and a related 
fixed point. Section \N\ discusses the spectrum of the 
Hamiltonian and explains how the cycle is related to a 
geometric sequence of binding energies (discrete levels) 
that approach zero. Then Section [V] discusses univer- 
sality and Wegner's marginal and leading irrelevant op- 
erators in a linearized approximation, with non-leading 
operators treated in Sect ion IvTl The RG analysis enables 
us then to show details in Section lVIII of how one can tune 
a Hamiltonian so that the limit cycle structure becomes 
rapidly visible in the spectrum. Section IVIII includes a 
numerical example with first two leading irrelevant oper- 
ators removed through the tuning so that the remaining 
irrelevant corrections to the limit cycle disappear at the 
rate given by factor 1/512 per cycle, for 6 = 2. Sec- 
tion describes generic properties of limit cycles and 
Section 17x1 briefly concludes the paper. 



The eigenvalue problem with H given by Eq. (JIJ, 



N 



(2) 



n— — oc 



can be solved for ip m , m < N, assuming that one knows 
E, by using Gaussian elimination. One solves for ipN in 
terms of all other components, ip n with n < N. Then, one 
expresses ipN-i in terms of components ipn with n < N — 
1, and so on. After the first p such steps, for negligibly 
small E, one is able to recognize the existence of limit 
cycles in the coupling constants g and h of period p or 
less, when they occur. 

The eigenstate components are re-written as ip n — 
6~"/ 2 </> ra for all n < N, and one defines 



N 

&N = ^ <Pn ) 

71= — QO 

N 
n— — oc 



(3) 
(4) 



The eigenvalue condition produces then a set of equa- 
tions, 



E 

E m 



9N&N - ihN^Nm — , (5) 



for all integers m < N. For m — N, ttnn = &n-i and, 
(1 - 9n - £a0 4>n = (gN + ihjf) cjv-i , (6) 

where = E/Ejy. Equation © gives 4>n in terms 
of all 4>m with m < N contained in ctjv-i- This result 
is then inserted into all the remaining equations with 
m < N — 1, leading to the new set of equations with the 
highest energy component removed, 

(1 - Cm) 4>m — gN&N-l ~ ihjfl^N-lm 

+ {~9n + ihtf) (gw + ih N ) (?n-i = , 



1 - 9N - ejv 



(1 - e m ) <j) m - ( g N + - 9n + fejv ] it v _| 



(7) 



1 - gN - ejv . 
- ih N ir N ^ lm = 0. (8) 



These equations appear to represent a new eigenvalue 
problem for a Hamiltonian matrix with elements, 

Hmn(gN~l, hN-l) = (EmEn) 1 ^ 2 

x (Smn — 9N-1 ~ ^JV-l s mii) , 

(9) 
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where 



9N-1 
flN-l 



9n - 

h N 
h. 



9% 



9n — 



(10) 

(11) 
(12) 



The new cutoff is Ajy-i = b N ^ = A/6. When N is 
reduced from N — 1 to A^ — 2, ejy goes over to ejv-i, and 
so on, along with the cutoff being reduced by successive 
powers of b. The coupling constant g^- P obtained after 
p such steps depends on gjv_ p +i and e^-p+i- 

Besides the sequence of coupling constants, the RG 
recursion formula can be used to find the eigenvalues E. 
A procedure to do so can be used to show the connection 
between the universal RG limit-cycle behavior and the 
energy spectrum. In principle, one specifies an energy E 
and computes the sequence of couplings g n for that value 
of E from the recursion 



-1 = 9r, 



1 - g n - E/E„ 



(13) 



starting from n — N and iterating all the way down to 
n = M+l, whereby one arrives at a gM that must satisfy 
the condition 



1 - g M - E/E_ 



A I 



0. 



(14) 



if E is an eigenvalue. Although one does not know in 
advance what E to pick, there exist simple search rou- 
tines to find out values of E that make the expression 
1 — gM — E/Em equal to zero. The solution to the re- 
cursion formula of Eq. (|13|1 will be developed in the next 
section. We will also employ later the reverse relation to 
Eq. (H3J), 



9n+l 



9n 



+ h 2 



Ej E n -\ 



(15) 



III. LIMIT CYCLE AND RESULTING FIXED 
POINTS 



Equation (|13fl can be rewritten as, 

g n ~i (g n /h) + h/(l - e n ) 



h " l-(g n /h)[h/{l-e n )} ' 

to exhibit the structure of a trigonometric identity, 

, . tana + tan/3 

tan (a + p) — . 

1 — tan a tan p 

Therefore, one can introduce the angles a n , and (3 n , 

9n 



a n — arctan 



(3 n = arctan 



h ' 
h 

l-£r. 



(16) 

(17) 

(18) 
(19) 



and obtain a simplified equation, 

a n _i = a n +fi n . (20) 
After p steps, one obtains 

ot n - v = a n + j(E, n, n - p + 1) , (21) 

where 

j(E, k,l) = p k + .... + A . (22) 

It is clear that the coupling constant returns to its value 
after the p steps, if 



j(E, n, n — p + 1) = 7r . 



(23) 



Although this is unlikely for arbitrary n and given E, 
a significant simplification occurs in the recursion when 
E — > 0, or, equivalently, when e n is so small that it can 
be neglected. In these circumstances, 



0„ 



arctan h 



1 + h 2 b n (1 + h 2 ) 2 



h(l-h 2 /3) (E 



(1 + h 2 ) 3 



(24) 



where 0(E 4 ) denotes terms order E 4 and higher. The 
simplification occurs because all angles (3 n become equal 
in the limit E — > 0, 



and 



(3 n = arctan h 

= p, 



/itan (a n + p/3) 

9n , 



(25) 
(26) 



(27) 



if P — n/p. Thus, one obtains a cyclic behavior of g n 
with period p , if 



h = tan (tt/p) 



(28) 



with an integer p, and if there exist eigenvalues equivalent 
to zero. The latter condition corresponds to the existence 
of two-body bound state at threshold in the three-body 
dynamics that was mentioned in the introduction. 

Equation l|28ll for h shows how the cycle emerges for 
integer p > 3 [2|. p = 2 requires infinite h, and p = 1 is 
equivalent to p = oo and h = 0. In that case, the cycle 
is infinite and one has asymptotic freedom instead of a 
finite cycle. For rational values of p in Eq. I|28l) . say, 
p = Pi jp2 with pi and P2 both integers, the sequence of 
coupling constants goes over a number of twists over the 
full period p^ > 3. For irrational values of p in Eq. 128fl . 
chaotic behavior appears. We limit further discussion to 
integer p > 3 in Eq. I|28|) . The sequence g n , g n -i, 
9Vi-p+ii repeats itself for as long as E is equivalent to 0. 
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Equation i|27|l has an important implication that the 
recursion formula must simplify when one computes g n - p 
directly from g n , because if E = 0, g n -p must equal 
g n regardless of the value of g n . This result suggests 
that one should study the relationship between functions 
fn(en) = g n (e n )/h and /„„ p (e„_ p ) = g n -p(e n -p)/h when 
E 7^ 0. Note that e„_ p = re„, where 



b". 



(29) 



is the energy scale-factor associated with moving down 
in energy from E n to E n _ p , which means over the entire 
cycle. But in order to derive the recursion that connects 
fn{x) with / n _ p (rcc), one needs to introduce p — 1 inter- 
mediate functions. 

Namely, we introduce functions f n -k{x) such that 
fn-ki^n-k) = g n -k(^n-k)/h, with e n - k = b k e n . In this 
notation, the recursion formula l|16|l reads, 



/n-l(^n-l) 



fn{x n ) + Zi(x n ) 

zi{x n ) 



(30) 



where x n -\ — bx n and z\{x) = h/(X — x). The RG 
transformation over two steps gives 



where 



z 2 {x) 



fn{x n ) + Z 2 (X„) 
1 - fn(x n ) Z 2 {x n ) 

Z\(x) + zi(bx) 
1 — z\{x)zi(bx) 



(31) 



Proceeding by induction, suppose that after k steps one 
has 



fn—k ip^n—k) 



1 fnip^ri) %k{%"n) 



(33) 



and one performs one more step according to Eq. 
The result is, 

it \ fn—k{,Xn — k)~\~Z\(x n — k) / <s 

Jn-fc-l^n-fc-l) = r — r . (34) 

1 - Jn-k{Xn-k) Zi{Xn-k) 

We substitute here the assumed expression (13*51 for 
fn-k(xn-k), and obtain, 

f ( \ /»(^») +^fc+l(^n) 

/n-fc-i(.^7i-fc-lj = i j ^ ^ 777^7 , (35) 



where 



1 - fn(x n ) Z k+ i(x n ) ' 

z k (x) + z^x) 
1 - z k (x) zi(b k x) 



(36) 



This sequence of intermediate steps implies that after 
p of them, 



j / \ /n(-^n) "I - (^n) 

1 "~ fnip^n) %p i%n) 



(37) 



where z p (x) is a function that can be calculated starting 
with Zi(a;) = h/(l — x) and using recursion given in Eq. 
(|36|l . The point is that one must have 



for the identity 



,(0) = 



/n-p(0) = /n(0) 



(38) 



(39) 



to hold with arbitrary values of f n (Q), which was ob- 
served before as a condition for the cycle to close in the 
case E = after p steps, Eq. (J27J. This implies that the 
Taylor series expansion for z p (x) around x — 0, starts 
with a term linear in x. 

To give an example, we take p = 3, the shortest pos- 
sible period. We focus on this example because when 
we discuss tuning to criticality for matrices smaller than 
100 x 100, the shortest period will allow us to display 
approach to the limit cycle spectrum most clearly. We 
have, 



fn(x) + Z 3 (x) 
1 - fn(x) Z 3 (x) ' 



where, 



Z3(x) 



a + ai + a 2 — a a 1 a 2 
1 — (a cii + a\a 2 + a 2 a ) 



(40) 



(41) 



(32) and a,j = h/(l — V x). For b — 2, one obtains 



Z3(x) 



(7/4)x(l - x) 



h 1 + x 3 - (7/4)x(l + x) ' 
which has an expansion, 

= (7/4)a; + (21/16).x 2 + (343/64).t 3 + 0(x 4 ) , 



(42) 



(43) 



that starts with a term order x, as expected. For x = 
0, i.e., for E = 0, one obtains a fixed point g n -3k = 
9n = 9n l° r arbitrary integer values of k. There are three 
such fixed-modulo-cycle points in the cycle of period 3, 

g n -3k = 9ni 9n-l-3k = and 3„_2-3fc = g£_ 2 - Tlic 

value of n does not matter. 



IV. THE BOUND-STATE SPECTRUM FOR 
THE LIMIT CYCLE 

In this section, we show that the limit cycle leads to a 
spectrum in the form of multiple exact geometric series 
with eigenvalues in each series separated by the factor r 
of Eq. Jill) in the limit N -> +oo and M -> -oo. One 
of these series describes bound states and is negative. We 
prove this by showing that for fixed E, M near enough to 
— oo, and N near enough to +oo, the sequence of coupling 
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constants g m (E) approaches, for m very negative but still 
well above M, 



g m {E)=-E{l-l/b)/b" 



(44) 



The above result can be obtained from Eq. (|15f) as- 
suming the eigenvalue condition qm{E) = 1 — E/Em, 
which means 



g M {E) = l-E/b 



M 



(45) 



The second term is huge in comparison to 1 or ft 2 for 
finite E and M — » — oo. Then Eq. I|15ll gives, 



g M +i(E) = g M (E) 



9li(E) 



l+g M (E)-E/b M + 1 



(46) 



By inspection, gM+i(E) = gM{E)/{\ + b) plus correc- 
tions of order 1 or ft 2 . If gM+ n -i{E) = —gM(E)/a n -i, 
then, neglecting terms order 1 and ft 2 again, gM+n{E) = 
—gM(E)/(a n -i + b n ), which implies that a n = 1 + b + 
... + b n and for finite m = M + n one obtains 



g M+n {E) = -E{l-l/b)/b M+ '> 



(47) 



as promised. But the existence of the limit cycle means 
that in the limit of large positive n the coupling constants 
return to the same value after p steps. 

Let gN — g* (E) in the limit as N — > oo through steps 
starting with a given finite m = M+n with large positive 
n and M — > — oo, and then sending k to infinity through 
steps m + p , m + 2p , m + 3p , . . . , m + kp , . . . , etc, thus 
avoiding other values on the cycle. Now suppose that E is 
increased by exactly r — b p . Equation (|47|l says that with 
k = 1 one has g m +p ( r E) — g m (E). But by translational 
invariance in m by p , which is valid in the limit cycle, this 
equation continues to hold for all larger values k. In the 
limit of large positive k, g m +k P (rE) = g m +(k-i) P (E) = 
g*(E). This means rE is also an eigenvalue for the same 
g* as E is. The same argument works for all values of m 
within one cycle, leading to as many series as elements 
in the cycle, which equals here p . 

It is necessary to actually construct p eigenvalues not 
separated by the factor r using Newton's method to de- 
termine the full set of v distinct geometric series. Our 
numerical experience |2| is that for any combination 
of b and p that we have tried, there is precisely one 
negative eigenvalue sequence and p — 1 positive eigen- 
value sequences. Besides the numerical experience with 
1 < b < 10, one can also consider very large b. In 
this case the eigenvalues are E n (l — g n ), which trans- 
lates, using Eq. (|18|l . into E n (l — fttana„), so that one 
might expect more negative eigenvalue series, one for ev- 
ery g n > 1 in one cycle [2| . But we observe here that with 
ft = tan/? a coupling constant g n can be greater than 1 
when (tt/2 — (3) < a n < 7r/2, and this can only occur for 
one n value in a cycle. This implies that there is only 
one negative eigenvalue series for all choices of 6, g* , and 
p, i.e., including the large b limit. The apparent discrep- 
ancy between this result and Eq. I|47|l disappears when 
one observes that Eq. P7|) is valid for \E\ ^> Em+u- 



The above analysis is based on the RG transforma- 
tion and attempts to avoid specific features of the model, 
besides the actual form of the transformation and the 
fact that we could use it in reverse. But our model pro- 
vides also an alternative path which takes advantage of 
the model simplicity and can be used as a crosscheck. 
Namely, one can generate energy eigenvalues as follows. 
First of all, for exceedingly large N, and E fixed, and M 
close to — oo, there is an eigenvalue whenever 

&n + Pn + 0N-1 + ■ • • + Pm = tt/2 + mix , (48) 
with m arbitrary. This means that 

0jv = fttaja[7r/2-7(-E,iV,-oo)], (49) 

where 

j(E, N, M) = arctan ■ 



arctan 



1 - EjE N _ x 



arctan 



1 - E/E. 
h 

1 - Ej E M 



(50) 

The sum becomes infinite when M goes to — oo but the 
sum converges because the arc tangents go to zero pro- 
portional to E n as n becomes more negative. Now we 
note the identity 



j(E, N-p, -oo) = j(rE, N, -oo) . 



(51) 



This means that E and rE can both be eigenvalues for 
the same g^ provided that j(E, N, — oo) = ~f(E, N — 
p, — oo). This is certainly true if E = and p arctan h = 
7T. But we can go further, dealing with the case of large 
but finite N. Suppose rE§ is an eigenvalue for a given 
gN- Then we can find an eigenvlue E near Eq from the 
requirement that 



7(£, N, -oo) = 7T + -f(rE ,N, -oo) . 
This becomes, using expansion of Eq. I)24[l. 

j(E, N, — oo) = 7(-Eo, N — p, — oo) +p arctan ft, 

ft I Eq t Eq Eq 

En 



(52) 



1 + h 2 \En-p+i En- p 



+ 2 



where 



N-p 



n=M 



E n [(l-E Q /E n ) 2 + h 2 ] 



+(E-E )S, (53) 



(54) 



From this result, which is valid to leading order in E 
(replaced by E ) and in E — Eq, respectively, we obtain 
the equation for E — Eq, 



E-Eq 



-h 



E, 



o 



1 + ft 2 \En- p +i 

Eq , , Eq \ 

En~) 



+ —^+...+-^)/S. (55) 

&N-P+2 i^N J 

This gives an E — Eq proportional to Eq/ En- p +i, which 
is the behavior expected from the presence of the leading 
irrelevant operator. 
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V. MARGINAL AND IRRELEVANT 
OPERATORS 

In this section we present an analysis of the two- 
variable renormalization group, in which f n -kp{x) of Sec- 
tion [ffi] is represented by a linear approximation in x for 
any k, 

fn-kp {Xn-kp ) = 9n/h + f^l kp X n - k p , (56) 

neglecting terms of order x 2 l _ kp or higher. The subscript 
n that numbers distinct points in the cycle will be omit- 
ted to simplify notation wherever possible. 

The RG equations for g* and for a jump of p steps 
are, 

9* = 9* , (57) 
rfi 1 l P =fi 1) + (l+9* 2 /h 2 )zW, (58) 

where Zp 1 ^ is the coefficient of x n in the expansion of 
z p (x n ) in powers of x n . The first step in the analysis is 
to solve Eq. I|58|) f° r a Hamiltonian with gn given and 
fixed, which means that fffl is zero. For N — n large, 
and a multiple of p , /« is fixed at a value which ensures 
that fn-p is the same as fn\ namely 

/« = /«* (59) 
= (l + . 9 * 2 /^ 2 )4 1) /(^-l), (60) 

which is valid for n of the form N — kp with k large 
enough. The solution leading from the boundary condi- 
tion fjip = to the fixed point is, 

fgl* = z^(l+ 9 * 2 /h 2 ) 

x (r -1 + r -2 + ... + 7- -Cfc-1) ) . (61) 

Our discussion of the fixed point and operators near 
the fixed point uses the variable x without any reference 
to n or k, and the function fix) is assumed to have the 
generic form, 

fix) = fo +hx. (62) 
The fixed point condition is, 

f(rx) = f(x) + (l + f 2 /h 2 )z^x. (63) 
Or, written in terms of a transformation R, 

fix) = R[f] 

= fix/r) + (l + f 2 /h 2 )z^x/r. (64) 

This has already been solved to give a one parameter 
family of fixed points, 

/*(<?*,*) = / * Of) + Kig*)x. (65) 



where 

/o(<?*) = 9*/h, (66) 

/ron = f {l) *ig*). (e?) 

Operators are obtained by linearizing the transforma- 
tion R about a fixed point. We write, 

f(x) = r(g*,x)+df(g*,x), (68) 

and treat df as small compared to /. Then, 

R[f* + df] = f(x/r)+df(x/r)+zW 

x [l + f 2 /h 2 + i2f Q df /h 2 )x/r}. (69) 

Wegner's operators are extracted from an eigenvalue 
condition, 

L(df) = wdf(x), (70) 

where w is an eigenvalue and L is the linearized form of 
R. Assuming a linear function df, 

df{x) — c + ci x , (71) 

one obtains two exquations that must be satisfied simul- 
taneously, 

c = w c , (72) 

- ci + 24 X) 9* - co = w Cl . (73) 
r r 

The eigenvalue problem has two solutions, one corre- 
sponding to a marginal operator with wq = 1, and one 
corresponding to an irrelevant operator with wi = r~ l = 
l/b p , which is less than 1. The marginal operator has Co 
arbitrary (let it be 1) but it has a c\ value also, which 
is determined given Co- The irrelevant operator has Co 
equal to zero and an arbitrary value for c\. Critical ex- 
ponents, A; for Z-th operator, can be read from Wegner's 
eigenvalues |3j] using relation, 

wi = r Xl . (74) 

Clearly, Ao is for the marginal operator. This operator 
is just the derivative of the fixed point function with re- 
spect to g*. Ai is — 1 for the irrelevant operator. This 
result is obtained here generically for all choices of pa- 
rameters in the model. 

The link from operators to the behavior of eigenvalues 
is following. The marginal operator causes the entire 
eigenvalue spectrum to shift while preserving the fact 
that successive eigenvalues have a ratio of r as long as 
eigenvalues are far from either the ultraviolet cutoff b N 
or the infrared cutoff b M . But as one approaches cither 
the ultraviolet or infrared cutoff limits, in the energy, the 
eigenvalue ratios change in ways that can only be de- 
termined numerically, except for the case that g* is cy- 
cled continuously through all real values from g* to co, 
then from — oo up to g* again. This complete cycle can 
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be shown to be equivalent to the renormalization group 
transformation through p steps already discussed. 

The leading irrelevant operator, if present (as it is if 
one solves the Hamiltonian of this article for finite N), 
causes the ratios of eigenvalues to approach the number r 
as the eigenvalue itself decreases from its maximum near 
En at a rate given by W\ = 1 jr for each reduction of the 
eigenvalue by r. This phenomenon occurs through the 
same recursion features that we used in Section IIVI to 
show the existence of geometric sequences of eigenvalues, 
except that one has to look for corrections due to finite 
N. If the leading irrelevant operator is tuned out (see 
Section |VTT|) . the next operator causes the rate to become 
1/r 2 . It will be shown in the next section that Wegner's 
irrelevant operators of successive orders have eigenvalues 
given by successive powers of 1/r, which suggests that 
the elimination of the next-to-leading irrelevant operator 
should accelerate the rate to 1/r 3 . Numerical example of 
this behavior is shown in Section lYlII 



VI. HIGHER-ORDER CONTRIBUTIONS 

In order to illustrate numerically in Section IVIII how 
the tuning of a finite initial Hamiltonian to criticality 
can be achieved, we need to discuss the behavior of the 
higher-order terms in f(x), i.e., the terms beyond fa and 
fi introduced in Eq. (|62[1 . Our numerical example in Sec- 
tion |HH uses N = 17 and M = —51, and we choose the 
smallest possible cycle period, p = 3, for the largest pos- 
sible number of cycles to fit within the bounds of these M 
and N. Accordingly, we focus in this section on the case 
p = 3 and discuss the higher-order terms in functions 
fn-3k(x), f n -3k-i{x), and f n ~3k~2{x), with arbitrary n 
and h = tan 7r/3 = a/3- 

In principle, the choice of g^ is not important, but 
in order to have the limit cycle universality pattern de- 
scribed in terms of integers or ratios of small integers, it 
is useful to aim at g* = — 7 at some value of n < TV, see 
Table H| and Appendix |X| 



TABLE I: The one parameter family of fixed points from Eq. 
igSJ for p = 3 and g* = -7. 



n 


hfS 


hfl 


n 


-7 


13 


n- 1 


-1/2 


13/16 


n-2 


5/3 


13/9 



For some n, the fixed point Eq. (|4U[) can be written in 
the form, 

/(8a:) = f{x) + z 3 [1 + /(8a:) f(x)} , (75) 

which can be used to generate a power series solution for 
the function f{x). The constant term, which we denoted 
g*/h, is arbitrary, but given g*, the remaining terms are 



determined. We can write 

hf(x) = hf(x,g*) 

= .9* + c x {g*)x + c 2 (g*) x 2 + c 3 {g*)x 3 + ... , 

(76) 

and the solution for the case g* = —7 is, 

hf(x,-7) = -7+ 13a;- (65/3)x 2 + (221/6)a; 3 

+ 0(z 4 ). (77) 

This result is extended to arbitrary values of b, p > 3, 
and g* , in Appendix 151 

Wegner's operators are infinitesimal functions df{x, g*) 
which satisfy Eq. (|40|l when added to the fixed-point 
solution for f*(x,g*). Namely, 

f(x,g*)+df(x,g*) + z 3 
l-z 3 [f(x,g*) + df(x,g*)} J{ ^ 9 ) 

= wdf(8x,g*), (78) 

valid to first order in df, with w being Wegner's eigen- 
value. Each eigenoperator begins with a different power 
of x, and the leading power of x determines w, namely 
Wj is 1/8 3 for the operator that begins with the power 
x^ , where j can be 1, 2, 3, etc. Each operator then has 
coefficients for all higher powers of x, just as f(x,g*) 
does. Key details of the supporting argument are given 
in Appendix IH1 for arbitrary b, g* , and integers p > 3. 

Extracting the linear terms, the precise equation for 
determining higher-order terms in df is the following, cf. 
Eq. O, 

wdf(8x,g*) = df(x,g*) + z 3 (x)[df(x,g*)f(8x,g*) 
+ f(x,g*)wdf(8x,g*)}. (79) 

VII. TUNING TO A CYCLE 

We show below how to tune parameters so that a finite 
Hamiltonian can exibit universal features of a limit cycle 
in its spectrum. We consider the case of p — 3 and we 
quote numerical results for b = 2. 

Our numerical studies began with alteration of the ma- 
trix element Hn n and attempts to use this element to 
obtain a fixed point value for the leading marginal op- 
erator already in two iterations. The idea was to obtain 
the approach to scaling with the rate of 1/64 instead of 
only 1/8 per cycle. Then, alterations of four matrix ele- 
ments with largest subscripts were studied, and it was 
found that by changing only these 4 matrix elements 
in the Hamiltonian of Eq. QJ, one can & l so make the 
next-to-leading irrelevant operator vanish, which leads 
to extremely fast approach to the geometric scaling of 
the eigenvalues, clearly visible within a small number of 
cycles. 

We look at mappings from parameters in the initial 
Hamiltonian matrix with additional couplings in the ele- 
ments -ffjv n, Hjsr jv-i, Hjsr-iN, /fjyjv-i) to the coupling 
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qn-2 and its expansion for small E. Denote the relevant 
matrix elements as, 



H 



N N 



Ext , 

H N -i n = \/E N E N -i(u + iv) , 
Hn N-i = V 'E N E N -i(u - iv) , 
Hn-i,n-i = Ejsi-xz . 



(80) 
(81) 
(82) 
(83) 



Then, gN-2 can be written as the ratio of two determi- 
nants. The numerator is the determinant of a 3 x 3 matrix 
that can be written with overall energy factors En, etc., 
removed. Then, the 3x3 matrix is, 



t — x/b 2 u — iv —g — ih 
u + iv z — x/b —g — ih 
—g + ih —g + ih —g 



(84) 



where g is gN and x is E / En-2- The denominator is the 
determinant of the upper left 2x2 matrix within this 
3x3 matrix. The result is that 

gN-2 = g - (g 2 + h 2 ) 

2u- (t + z) + x(b+ l)/b 2 



tz - (u 2 + v 2 ) - x(tb + z)/b 2 + x 2 /b 3 ' 

(85) 

It is convenient to set t = u = z, leaving z and v as free 
parameters. Then, with h = 

gw-2 = g + x(g +3) —^r 



x 2 (g 2 + 3) z 



b 2 y 2 



+ 0{x 3 ). (86) 



The fixed point value for gN-2, through order x 2 , can 
also be worked out, and gives 



gN-2 = g + x (g 2 + 3) 



1 



4(6-1) 

+ * 2 tf + Z) 9+( '- m l +l) +0 { x*). 



16(6-1) 2 

From matching these two formulae, one obtains 
v 2 =4(6 2 - l)/6 2 = 3, 

and, 

z = -g- (6-l)/(6+l) = 20/3, 



(87) 
(88) 
(89) 



with the numerical values valid for g = —7 and 6 = 2. 

Using t = u = z = 20/3, v = ^/3, and gN — —7, we 
find that 



9N-2{x) 



-7+ 13x/(l + 5a;/3-a; 2 /24) 
65 

y 



-7 + 13a; 



■ x 2 +o(x 3 ) , (90) 



TABLE II: The left column contains eigenvalues of the initial 
Hamiltonian which is tuned so that the leading and next-to- 
leading irrelevant operators reach already at iV — 2 (see 
the text for details, t = u = z = 20/3, v — \/3, = —7, 
p = 3, b = 2, N = 17, M < -51). The eigenvalues are 
ordered by magnitude, and the negative ones are separated 
from the positive. The middle column contains the ratio to 
prior (smaller in magnitude) eigenvalue. The right column 



indicates the factor w.. 



r with which corrections to the 



ratio r — 8 disappear when binding energies approach 0. The 
number ~ 400 in the right column turns into 511.6 when 
M = -51 is changed to M = -60. 



E 


ratio 


w 3 


0.115359460521 






0.326274994387 


2.828333219599 




0.922875684170 


2.828521033011 




2.610199955097 


2.828333219599 




7.383005473358 


2.828521033011 




20.881599640790 


2.828333219600 




59.064043787554 


2.828521033043 




167.052797170395 


2.828333220314 




472.512353132001 


2.828521049247 




1336.422560688180 


2.828333591344 




3780.110931254790 


2.828529719903 




10692.223889183900 


2.828547649430 




30306.851247105530 


2.834475929536 




91327.456629399810 


3.013426102394 




1824274.370928830000 


19.975092247797 




-3.321031582839 






-26.568252662770 


8.000000000017 


~ 400 


-212.546021531960 


8.000000008649 


511 


-1700.369094987040 


8.000004341325 


502 


-13606.149375382800 


8.001879953885 


433 


-111890.939577163000 


8.223556605927 


119 



in agreement with Eq. I|77fl . and indeed the eigenvalues 
converge very fast towards a fixed ratio, as illustrated by 
Table 1TT1 There is a negative eigenvalue sequence sepa- 
rated by a factor rapidly approaching 8 and a positive 
eigenvalue sequence whose ratios rapidly approach an al- 
ternating sequence of two ratios both near the square root 
of 8. As a numerical check, we have verified that the sum 
of these eigenvalues is to twelve decimal places the trace 
of the original Hamiltonian. The eigenvalues in Table 
|n]are given for N = 17 and M = — 51 (low enough to 
have no impact on these numbers: smaller eigenvalues are 
omitted). The number ~ 400 in the right column turns 
into 511.6 when M = —51 is changed to M = —60. This 
result merely confirms that the numbers in the table cor- 
respond to the limit cycle with M — *■ — oo. The ratios for 
the positive eigenvalues can be compared with b p >( p ~ ^ = 
V8 - 2.82842712474619, looking at the average value 
(2.828333219599+2.828521033011)/2 = 2.828427126305, 
and at the product 2.828333219599 x 2.828521033011 ~ 
7.999999999999491. 
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VIII. GENERIC PROPERTIES OF LIMIT 
CYCLES 

The model Hamiltonian discussed in this paper has an 
extraordinarily simple solution despite its exhibiting a 
limit cycle. The initial coupling g^ is given as an explicit 
analytic expression as a function of an energy eigenvalue 
E that involves nothing worse than a convergent infinite 
sum of arc tangents. The renormalization group transfor- 
mation for the model can be written in terms of one, two, 
or an arbitrary number of coupling parameters. Analy- 
sis of the transformation with j couplings yields exact 
results for the limit cycle itself and for the critical ex- 
ponents of the j leading operators, including a marginal 
operator and a succession of irrelevant operators with in- 
teger exponents. 

Obviously, this is too simple to all carry over to more 
general renormalization group limit cycles. What does 
carry over to the more general case? First of all, we can- 
not guarantee that a renormalization group limit cycle 
will lead to any bound states. But if there are bound 
states close to threshold, then these bound states must 
form an infinite geometric series converging on zero. The 
reason for this infinite geometric series is the existence of 
a fixed point Hamiltonian for a change of scale r that cor- 
responds to a transformation around the complete limit 
cycle (r — b p in the discrete model). The proof of this 
statement begins with the fact that if the fixed point 
Hamiltonian with cutoff A has an energy eigenvalue E 
much less than A, then the Hamiltonian after the renor- 
malization group transformation has cutoff A/r, yet ex- 
actly the same form. Therefore the Hamiltonian after 
the transformation must have an eigenvalue E/r. But 
the definition of the transformation is that it preserves 
eigenvalues, which means that the Hamiltonian before 
the transformation must have E/r as an eigenvalue too, 
as well as its original eigenvalue E. Then by following 
the same reasoning over again, one can establish that 
the complete geometric series E, E/r, E/r 2 , etc., must 
all be eigenvalues for the fixed point Hamiltonian. 

The second observation to make is that marginal and 
irrelevant operators vary at different locations on the 
limit cycle, but the critical exponents of these operators 
do not vary. In the model, the operators have the form 
of a set of coefficients ci(g*) , c 2 (g*), etc., that depend 
on the values that make up the limit cycle. The crit- 
ical exponents are the integers —1, —2, etc., and do not 
depend on the location n on the cycle. Let us us then 
consider the general case. Suppose there is a non-linear 
transformation Rb that forms the basic (this is why we 
introduce the subscript b) renormalization group trans- 
formation that when iterated p times, Rb Rb ■ ■ ■ Rb, pro- 
duces the transformation R defined earlier that corre- 
sponds to a complete circuit of the cycle, which we could 
denote here by R r . Suppose then that one has a limit cy- 
cle of Hamiltonians H* of period r — b p . Suppose that 
dH n is one of Wegner's operators for this limit cycle. We 
now apply the renormalization group transformation R r 



to H* + dH n k times, where k is arbitrarily large. Sup- 
pose that the critical exponent for dH n were an exponent 
A„ dependent on n, i.e., we would have different expo- 
nents for different ns. Then there would exist two ways 
to apply the iterated transformation that would give two 
different results, while on the other hand they would have 
to give the same result in the case of limit cycle. The first 
way is to replace Rb iterated kp times by R r iterated k 
times, in which case H* + dH n becomes H* + r Xn dH n . 
The second way is to apply Rb to H* + dH n , then ap- 
ply the transformation R for k — 1 times, then apply Rb 
for p — 1 times. The only source of k dependence is the 
k—1 applications of R r , but it is applied to H*_ 1 +dH n ^i 
which means that dH is multiplied by 7-( fc-1 )*n-i , if A„_i 
is different from X n , there is an immediate contradiction 
in the fc-dependence of the two different approaches. The 
same argument works for all values of n within the cycle. 

The third observation is that the model has a marginal 
operator, that is, an operator with a critical exponent of 
zero. The presence of a marginal operator is guaran- 
teed if one has a renormalization group in differential 
form so that the limit cycle must be continuous, and 
position on the cycle is labeled by a continuous param- 
eter, just as g* is continuous in the model. In this case, 
the marginal operator is obtained by differentiating the 
limit cycle Hamiltonians themselves with respect to the 
parameter , and it is easily established that the criti- 
cal exponent for this operator must be zero. But when 
the renormalization group is discrete, and the limit cycle 
is a finite set of points rather than a continuous curve 
in the space of Hamiltonians, we do not know whether 
or not a marginal operator has to occur. In addition, 
while our model has no relevant operators, this need not 
be true for other limit cycles, and indeed, the limit cy- 
cle for the three-body case has a relevant operator asso- 
ciated with infinitesimal changes in the strength of the 
two-body potential in the model. When the relevant op- 
erator is absent, there is a two-body bound state exactly 
at threshold, and inclusion of the relevant operator causes 
this threshold bound state either to acquire a non-zero 
binding energy or to blend into the two-body scattering 
continuum. 

Finally, we showed that parameters in an initial Hamil- 
tonian could be tuned to eliminate the two leading irrele- 
vant operators and thereby ensure a very rapid approach 
to a pure geometric series of the bound state eigenval- 
ues for the Hamiltonian. We cannot guarantee that such 
tuning will be possible in other examples of renormaliza- 
tion group limit cycles, such as the limit cycles for the 
three-body Hamiltonians 0, 0, 0, IS EH • The problem is 
that there is a non-linear mapping from parameters in an 
initial Hamiltonian to coefficients of Wegner's operators, 
and with non- linear mappings, there is no assurance that 
any given set of output parameters can be obtained from 
the mapping. We were lucky to find an easily identified 
solution for the tuning problem for the starting Hamil- 
tonian that we used. In other cases, one has to make 
careful studies to establish if a desired option for tuning 
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can work. 



IX. 



CONCLUSION 



The renormalization group limit cycle discussed in this 
paper is extraordinarily simple to analyze thanks to the 
analytic solution to the recursion formula for g n and to 
the analytic formula expressing cjm in terms of any eigen- 
value E. It is far simpler to study than the three-body 
Hamiltonian of references 0, IE IE IS [JTJl - which more 
than compensates for its lack of immediate practical ap- 
plications. We do not know how many other examples of 
Hamiltonians with renormalization group limit cycles will 
be found in the future. But hopefully, the clarification 
of the expected behavior of renormalization group limit 
cycle behavior described in this article will make it easier 
to search for and make sense of limit cycle behavior in 
more complex circumstances. Moreover, the possibility 
of bound states forming geometric series converging on 
zero energy provides a startling characteristic to search 
for. 



APPENDIX A: CALCULATION OF ENTRIES 
FOR TABLE H FOR INTEGERS p > 3 



Suppose that, 



9n = 9 n + c„e„ , 



where g* — hf^, and c„ = hf\. We have, 

ffn-l * +Cri-ie n _l = <?* + c„e„ 

{g* n + c„e„) 2 + h 2 



(Al) 



(A2) 



Expanding the right-hand side in powers of e n , and taking 
into account the cycle relation (|TT)|> between the coupling- 
constants g^-\ and <?* for E = 0, the result is, 



with 



Cn — 1 3<n ' Vn Cn j 



1 



Vn 



1 l + h 2 

b (T 



-,*\2 



(A3) 

(A4) 
(A5) 



<7* runs over a cycle and h and b are constants. Therefore, 
c n also approaches a certain cycle when n is lowered. This 
can be seen after p steps, 



with 



C-n—p T" V n C n , 

^n 3<n— p +1 Vn — p +1 ^n— p +2 

+ Vn-p+1 Un-p+2 x n-p+3 + • 

+ Vn-p+1 Un-p+2 ■ ■ ■ Vn-l %n i 

V n = Vn-p+1 Vn-p+2 ■ ■ ■ Vn-l Vn ■ 



(A6) 



(A7) 
(A8) 



The fixed point of this transformation over the cycle is 
given by 



u n /(l - v n ) ■ 



(A9) 



and the p — 1 values in between can be calculated for 
given g n at one value of n using Eq. i|A3|l . 



APPENDIX B: DIRECT ENERGY EXPANSION 
IN THE MODEL 

Including corrections due to E ^ in Eq. I|24|) . and 
inserting them in Eq. I|21() , one obtains for p arctan h = it 
that, 



^-^ = tan [a* + e x E + e 2 E 2 + e 3 E 3 + 0(E 4 ) 



where g n and a n are functions of E, and 



('i 



('2 



e3 



1 + h 2 b n ' 

h b 2 
(1 + h 2 ) 2 b 2n ' 

h ( _ hT 
(1 + h 2 f [} ~ 3 J ( V ; « 



(Bl) 

(B2) 
(B3) 

(B4) 



with bk — {b kp — l)/{b k — 1). Then, the coupling constant 
dependence on E obeys the recursion, 



9n—p 

ei 
e 2 



g n + e x E + e 2 E 2 + e 3 E 3 + 0{E A ) , (B5) 
g 2 n + h 2 h 



l + h 2 b n 
9l + h 2 

(i + h 2 y 



i b2 + -9» b i] -jpn- - 



e3 



9n + h 
(l + h 2 ) 3 



1 - y ) b 3 + 2g n b 2 b 1 



+ (9l + h 2 mbt] ^ 



(B6) 



(B7) 



(B8) 



Assuming that f n (rx) = g n - P {rE) and g n (E) = f n (x) 
with x — E/E n , and expanding Eq. 1|B5|) in a series of 
powers of x, one arrives at (the subscript n is dropped, 
but there are p distinct functions with p different values 
of g* , forming a cycle) , 

f{x) = g* + c[ x + c* 2 x 2 + 4 x 3 + 0(x 4 ) , (B9) 

with the coefficients, 

* g* 2 + h 2 h 



l + h 2 r - 1 

*2^ h 2 



[l + h 2 f 



9 



bi 
r-1 



(BIO) 
, (Bll) 
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g* 2 + h 2 
(l + h 2 f 

h 



2.9 



(g* 2 + h 2 /3) 

l>2 



1 r 2 



1 



r - 1 

(1 - h 2 /3)b. 
r 3 - 1 



(B12) 



These expressions provide fixed-modulo-cycle values of 3 
additional coupling constants for arbitrary choices of g*, 
h = tan(7r/p), and b, with integer p > 3, in addition 
to the results obtained in Section The p different 
constants g* in the limit cycle, g*, g*_ l% . . ., 
are determined by a choice of any one of them with an 
arbitrarily selected n. 



have fixed point values given by T m (g*) — (<?*), with 



T* m {9*) = l m /{r m 1) , 



(C5) 



where 7 m are given by the expansion of the known func- 
tion j(x), 



j(x) — it + ^2 lmX m /m\ 



(C6) 



rn — 1 



APPENDIX C: EIGENVALUES Wj FOR j > 1 

Equation (|21|l implies 

On- p (rx) = a„(i) + 7(1) , (CI) 

where 

p-i h 
7 (x) = arctan — y- , (C2) 

fe=0 

see Eq. (|22|) . One has 7(0) = n, and the fixed point 
relation can be written as, 

a* (rx, g* ) = a* (x, g*) + 7(1) - 7 (0) . (C3) 

The coefficients of the Taylor series for the function 
a(x,g*) around zero, 

00 

a(x,g*) = r (g*) + ^ r m {g*)x m /m\ , (C4) 

m— 1 



Suppose that 

a(x, g*) = a* (x, g*) + SjX J , (C7) 

with infinitesimal Sj and j > 1. The RG transformation 
over the cycle gives, 

= '' ''^ , (C8) 

and a whole series dfj (x,g*) of powers of x greater than or 
equal to j is generated, by Wegner's eigenvalue condition, 

Wjdfj(rx,g*) = dfj(x,g*)+z p (x)[dfj(x,g*)f(rx,g*) 
+ /!.<•.'/•:«', '//>••<••//•: • (C9) 

with f(x,g*) = h tan [a* (x, g*)]. 
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